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ABSTRACT 

d We use a free-energy minimization approach to describe in simple and clear phys- 

^ ical terms the secular and dynamical instabilities as well as the bifurcations along 

equilibrium sequences of rotating, self-gravitating fluid systems. Our approach is fully 
£S) nonlinear and stems from the Landau-Ginzburg theory of phase transitions. In the fi- 

nal paper of this series, we examine higher than second-harmonic disturbances applied 
to Maclaurin spheroids, the corresponding bifurcating sequences, and their relation to 
^ nonlinear fission processes. 

I/*-) The triangle and ammonite sequences bifurcate from the two third-harmonic neu- 

tral points on the Maclaurin sequence while the square and one-ring sequences bi- 
furcate from two of the three known fourth-harmonic neutral points. The one-ring 
sequence has been analyzed in Paper II. In the other three cases, secular instability 
Q | does not set in at the corresponding bifurcation points because the sequences stand 

^ and terminate at higher energies relative to the Maclaurin sequence. Consequently, 

;_i an anticipated (numerically unresolved) third-order phase transition at the ammonite 

c/2 bifurcation and numerically resolved second-order phase transitions at the triangle and 
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square bifurcations are strictly forbidden. Furthermore, the ammonite sequence exists 
at higher rotation frequencies as well and is similar in every respect to the pear-shaped 
sequence that has been analyzed in Paper III. 

There is no known bifurcating sequence at the point of third-harmonic dynamical 
instability. This point represents a discontinuous A-transition of type 3 that brings a 
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Maclaurin spheroid on a dynamical time scale directly to the binary sequence while the 
original symmetry and topology are broken in series. The remaining fourth-harmonic 
neutral point also appears to be related to a type 3 A-transition which however takes 
place from the lower turning point of the one-ring sequence toward the starting point 
and then on toward the stable branch of the three-fluid-body (triple) sequence. A third 
type 3 A-transition, taking place from the one-ring sequence toward the starting point 
and then on toward the stable branch of the four-fluid-body (quadruple) sequence, is 
also discussed. 

The two-ring sequence bifurcates from the axisymmetric sixth-harmonic neutral 
point on the Maclaurin sequence also toward higher energies initially but eventually 
turns around and proceeds to lower energies relative to the Maclaurin sequence. The 
point where the two sequences have equal energies represents a fourth type of A 
transition which is not preceded by a first-order phase transition. This type 4 A 
transition results in double fission on a secular time scale: a Maclaurin spheroid breaks 
into two coaxial axisymmetric tori that rotate uniformly and with the same frequency. 

Finally, our nonlinear approach easily identifies resonances between the Maclaurin 
sequence and various multi-fluid-body sequences that cannot be detected by linear 
stability analyses. Resonances appear as first-order phase transitions at points where 
the energies of the two sequences are nearly equal but the lower energy state belongs 
to one of the multi-fluid-body sequences. Three nonlinear resonances leading to the 
turning points of the binary, triple, and quadruple sequences are described. 

Subject headings: galaxies: evolution - galaxies: structure - hydrodynamics - insta- 
bilities - stars: formation 

1 Introduction 

In this series of papers, we explore the physics of bifurcation and instability points 
along equilibrium sequences of uniformly rotating, self-gravitating, incompressible fluids 
using a method that stems from the Landau-Ginzburg theory of phase transitions (see 
Landau & Lifshitz 1986; see also related works by Bertin & Radicati 1976, Lebovitz 1977, 
Poston & Stewart 1978, Constantinescu, Michel, & Radicati 1979, Gilmore 1981, Hachisu 
& Eriguchi 1983, 1984a, Eriguchi & Hachisu 1983b, 1984, Tohline 1985, and Christodoulou 
& Tohline 1990). In the first three papers (hereafter referred to as Papers I III) , we have 
described the differences between the second-order phase transitions that correspond to the 
second-harmonic instabilities on the Maclaurin sequence (Paper I), we have investigated the 
first-order phase transitions and the discontinuous A-transitions related to the one-ring, 
dumbbell-binary, and Maclaurin toroid sequences (Paper II), and we have analyzed the 
third-harmonic bifurcation of the pear-shaped sequence from the Jacobi sequence and its 
relation to the classical fission hypothesis (Paper III). 

Phase-transition theory allows us to understand all points of bifurcation and instability 
along equilibrium sequences as points where phase transitions of various orders may appear 
and may either be allowed or may be forbidden depending on the (non)-conservation of 
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four "integrals of motion." These integrals of motion are the mass, the energy, the angular 
momentum, and the circulation. 

The most important contribution of this series of papers is the realization that the above 
integrals of motion determine alone whether/where points of secular and dynamical insta- 
bility exist. In our terminology, integrals of motion determine whether a phase transition 
appears at a point and whether a known phase transition is allowed or forbidden. This means 
that an accurate description of an evolutionary path requires knowledge of whether any in- 
tegral of motion varies during evolution and understanding of how such an integral may vary 
in time. Consider two examples that are particularly relevant to the results described below. 
(1) In general, circulation is not conserved between sequences of uniformly rotating equilibria 
and poses no constraint to the occurence of the corresponding phase transitions. Even so, 
not every neutral point is a point of secular instability because a bifurcating "daughter" se- 
quence may unfold, at least initially, toward higher energies relative to its "mother" sequence 
(cf. Papers II and III). (2) The instabilities/transitions detected by linear stability analyses 
are not the only allowed phase transitions between sequences. In particular, linear analyses 
do not capture nonlinear resonances between sequences (see below), first-order phase tran- 
sitions that represent nonlocal changes, or type 2 A-transitions that represent discontinuous 
catastrophes (see Paper II). 

In this paper, we complete our analysis of the known critical points and the correspond- 
ing bifurcating sequences of the Maclaurin sequence. Specifically, we describe the phase 
transitions that are related to: (a) the two third-harmonic neutral points and the corre- 
sponding bifurcating (triangle and ammonite) sequences; (b) the third-harmonic dynamical 
instability point; (c) the three fourth-harmonic neutral points and the two known bifur- 
cating (square and one-ring) sequences; (d) the axisymmetric sixth-harmonic bifurcation 
of the two-ring sequence; and (e) nonlinear resonances between the Maclaurin and various 
multi-fluid-body (binary, triple, quadruple) sequences. 

In what follows, we consider for simplicity only the evolution of incompressible, uni- 
formly rotating fluid-masses that collapse under the action of self-gravity. We also adopt 
again the assumptions made in Papers II and III: (1) Mass is always strictly conserved not 
only between the initial state and the final equilibrium state but also in the intermediate 
states that an evolving fluid-mass passes through. (2) Energy always decreases as a fluid- 
mass evolves toward a new equilibrium state. (3) Angular momentum, expressed in physical 
units, is always strictly conserved but, in normalized units, it increases or stays constant in 
time. This general increase in normalized angular momentum j reflects a general increase 
of the density p during the evolution of a contracting fluid-mass [j oc p 1 / 6 in equation (1.2) 
below]. (4) Circulation is strictly conserved only in the complete absence of "viscosity" (i.e. 
in a "perfect" fluid) but its absolute value decreases during the evolution of a "viscous" 
fluid-mass. 

Our interpretations of instability points rely on numerical results obtained by Eriguchi, 
Hachisu, and Sugimoto between 1981 and 1986. We refer below to particular results from 
these numerical computations as we need them. We also adopt for convenience the dimen- 
sionless quantities u (rotation frequency), j (angular momentum), and E (total energy) that 
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were used by these authors as well as in Papers II and III. These quantities are defined by 

o 2 

,2 _ m2f " ^ 



^ 10 te). (»■») 



and 

where is the rotation frequency, G is the gravitational constant, p is the mass density, L is 
the angular momentum, and M is the mass, all in physical units. The term T + W denotes 
the total energy in physical units as the sum of the total kinetic energy due to rotation T 
and the total gravitational potential energy W. As we have seen in Paper I, the total energy 
T + W and its dimensionless counterpart E also represent the free-energy function when they 
are not constrained by equilibrium conditions (for more details see Tohline & Christodoulou 
1988). Finally, when u, j, and E refer to equilibrium configurations specifically we write 
them as u , j , and E , respectively. 

The remainder of the paper is organized as follows. In §2, we describe the triangle 
and ammonite sequences that bifurcate from the two third-harmonic neutral points on the 
Maclaurin sequence and the square sequence that bifurcates from one of the fourth-harmonic 
neutral points. These three neutral points are not points of secular instability. The second- 
order phase transitions toward the triangle and square sequences are forbidden because these 
sequences exist only at higher energies relative to the Maclaurin sequence. The ammonite 
sequence is also unstable and is further similar to the pear-shaped sequence that was an- 
alyzed in Paper III. All three sequences terminate prematurely because of equatorial mass 
shedding. Consequently, there exist no related first-order phase transitions or A-transitions 
similar to those found in Paper II to be associated with the one-ring sequence. In §3, we 
describe the two-ring sequence that bifurcates from the axisymmetric sixth-harmonic neu- 
tral point on the Maclaurin sequence. This sequence initially unfolds toward higher energies 
but eventually crosses below the Maclaurin sequence. A new type of symmetry-preserving, 
topology-breaking A-transition appears at the point where the two sequences have equal en- 
ergies. Because of the particular structure of the two-ring sequence, this type 4 A-transition 
is not preceded by a first-order phase transition. 

In §4, we describe the A-transitions of type 3 that appear to be associated with the 
third-harmonic dynamical instability point and with two of the fourth-harmonic neutral 
points on the Maclaurin sequence. The latter two transitions originate from points on the 
lower energy toroidal branch of the one-ring sequence. Type 3 A-transitions lead to nonlinear 
fission of a Maclaurin spheroid or of an axisymmetric torus to multiple fluid bodies while first 
the symmetry and then the topology of the original figure are broken. The initial breaking of 
the symmetry (which is characteristic of only type 3 A-transitions) must be the reason that 
this type of dynamical instability is detected by linear stability analyses (see Chandrasekhar 
1969, hereafter referred to as EFE). In §5, we describe the nonlinear resonances that exist 
between the Maclaurin sequence and various multi-fluid-body sequences. Resonances are 
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FIGURE 1. The rotation frequency squared (left panels) and the free energy E (right 
panels) along the triangle, square, ammonite, and Maclaurin sequences are plotted as func- 
tions of the dimensionless angular momentun squared j%. The Maclaurin sequence and the 
exact locations of the bifurcation points B (solid circles) have been calculated analytically 
(see also Table 1). All points (open circles) of the bifurcating sequences were obtained from 
Tables II IV of Eriguchi & Hachisu (1982) who gave the points denoted by the letter T as 
the termination points of the sequences due to equatorial mass shedding. 
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related to the turning points of the binary, triple, and quadruple sequences and appear 
as first-order phase transitions that also result in nonlinear fission for sufficiently strong 
disturbances. In §6, we summarize our results and we relate them to previous linear stability 
analyses and numerical hydrodynamical simulations. 

2 The Triangle, Square, and Ammonite Sequences 

The triangle, square, and ammonite sequences bifurcate from the Maclaurin sequence 
at the neutral points that correspond to the harmonic indices (/,m) = (3,3), (4,4), and (3,1), 
respectively (see EFE and Hachisu & Eriguchi 1984c). The physical properties of these 
sequences were computed by Eriguchi & Hachisu (1982). Their results are shown in Figure 
1 where the letters B denote the exact locations of the bifurcation points (EFE) and the 
letters T denote the termination points of the sequences due to mass shedding. In the E ) 
plane, the triangle sequence terminates above but close to the Maclaurin sequence while the 
square sequence turns around toward lower j -values and slightly higher .E^, -values before 
terminating. 

We see then from Figure 1 that there is no instability of any kind related to the triangle 
and square sequences. Consequently, the phase transitions at the bifurcation points of these 
sequences are strictly forbidden. Notice, however, that the slopes dE /dj appear to be 
continuous while the slopes duj /dj have a finite discontinuity at the bifurcation points. 
(This conclusion is certainly valid despite the mismatch seen at the bifurcation points in 
Figure 1 between the analytical and the numerical results.) This is consistent with our 
description of second-order phase transitions given in Paper I. 

The ammonite sequence is different than both the above sequences and similar in every 
detail to the pear-shaped sequence that bifurcates from the Jacobi sequence (Paper III). 
In this well, the sequence exists only at higher rotation frequencies and at higher 

energies relative to the Maclaurin sequence and no instability of any kind appears. We 
believe then that the ammonite sequence bifurcates in the (j%,u}g) plane smoothly toward 
lower j -values and higher u; -values initially and that the bifurcation point represents a 
forbidden third-order phase transition (compare the relevant panels in Figure 1 to the two 
figures in Paper III). This conclusion about the overall structure of the ammonite sequence is 
identical to the results of Paper III concerning the (3,1) neutral point on the Jacobi sequence 
and the behavior of the corresponding bifurcating pear-shaped sequence. 

For convenience, we list in Table 1 the physical properties of the above discussed neutral 
points as well as of the critical points discussed below. The meridional eccentricity e of the 
corresponding Maclaurin spheroid at each point is also listed along with the values of 
uj 2 01 E , and the ratio T/|W|. The second-harmonic instability points discussed in Paper I 
are also listed in Table 1 for completeness. The first entries of the (4,0) neutral point were 
calculated using the value e=0. 98531 (Chandrasekhar 1968). The second entries of the (4,0) 
neutral point are taken from Paper II (they were obtained from the results of Chandrasekhar 
1967) and are listed here for comparison purposes. The numerical differences between these 
two determinations are practically negligible. 
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Table 1 

Critical Points on the Maclaurin Sequence 



Point (I, m) 


Bifurcating Sequence e j% 




E 


T/\W 


SECOND-HARMONIC POINTS 


Neutral (2,2) 
Dynamical Instability^) 


Jacobi 0.81267 0.45547 
x = +1 Riemann 0.95289 1.27957 


9.35574 
11.00550 


-2.94820 
-6.36714 


0.13753 
0.27382 


THIRD-HARMONIC POINTS 


Neutral (3,3) 
Dynamical Instability^ 
Neutral (3,1) 


Triangle 0.89926 0.78526 
0.96696 1.53796 
Ammonite 0.96937 1.59504 


11.00366 
10.49141 
10.35347 


-4.55015 
-7.10294 
-7.24912 


0.20233 
0.30308 
0.30895 


FOURTH-HARMONIC POINTS 


Neutral (4,4) 
Neutral (4,2) 
Neutral (4,0) 


Square 0.93275 1.03712 
0.98097 1.96741 
One-ring 0.98531 2.17884 


11.23066 

9.33478 

8.71210 


-5.55125 
-8.07828 
-8.46563 


0.24174 
0.34290 
0.35925 


Neutral (4,0) 


One-ring (Paper II) 0.98523 2.17412 


8.72609 


-8.45757 


0.35891 


SIXTH-HARMONIC POINTS 


Neutral (6,0) 
Secular Instability^) 


Two-ring 0.99375 2.91500 
0.99802 3.98856 


6.61108 
4.21544 


-9.46087 
-10.27277 


0.40344 
0.44322 



(a) Also critical point of a second-order phase transition (Paper I). 

(b) Also a type 3 A-point (§4). 

(c) Also a type 4 A-point (§3). 



3 The Two-Ring Sequence 

The two-ring sequence is shown in Figures 2 and 3; it bifurcates from the Maclaurin se- 
quence at the (6,0) neutral point that corresponds to a meridional eccentricity of e=0. 99375 
(Hachisu & Eriguchi 1984c). The physical properties of this sequence were computed by 
Eriguchi & Hachisu (1982, 1983a). Like the one-ring sequence of Paper II, the two-ring 
sequence bifurcates smoothly in both Figures 2 and 3 and unfolds toward higher energies 
initially. This implies that the third-order phase transition at the bifurcation point is for- 
bidden as in the case of the one-ring sequence. 
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FIGURE 2. The rotation frequency squared along the two-ring and the Maclaurin se- 
quence is plotted as a function of the dimensionless angular momentun squared The 
Maclaurin sequence has been calculated analytically. The bifurcation point B (filled cir- 
cle) has also been calculated analytically using an angular momentum value of j% = 2.915 
(Hachisu & Eriguchi 1984c and Table 1). All points (open circles) of the two-ring sequence 
other than the bifurcation point were obtained from the corresponding tables of Eriguchi & 
Hachisu (1982, 1983a). The dotted line XC at f = 3.98856 denotes a type 4 X-transition 
(secular instability) that is not preceded by a first-order phase transition. 

In contrast to the one-ring sequence, four turning points can be seen in the two-ring 
sequence of Figure 2. These turning points are associated with changes in the structure of 
the equilibrium objects (Eriguchi & Hachisu 1983a) and cause a difference in the behavior of 
the sequence in the (j%,E ) plane of Figure 3: The two additional turning points make the 
sequence turn forward toward higher j -values. The two-ring sequence eventually crosses 
below the Maclaurin sequence but a first-order phase transition similar to that of the one- 
ring sequence in Paper II is not realized because there is no intermediate branch to act as 
a free-energy barrier at the point where the two sequences have equal free energies. Thus, 
this point is a degenerate critical point of the free-energy function. 
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FIGURE 3. The free energy E is plotted versus j% for the sequences of Figure 2. The 
intersection X=C between the Maclaurin and the two-ring sequence is a X-point of type 4 
(see also Table 1 ). 

The degenerate critical point (X=C) between the two sequences in Figure 3 occurs at 
j„ = 3.98856, E Q = —10.27277 and corresponds to point X on the Maclaurin sequence of 
Figure 2 with e = 0.99802, lo 2 q = 4.21544, and T/\W\ = 0.44322. This point is a A-point 
of a new type which we call type 4. The type 4 A-transition from point X to point C in 
Figure 2 results in axisymmetric nonlinear fission of the original Maclaurin spheroid to two 
toroidal figures that rotate uniformly and with the same frequency. Thus, this transition 
may only take place on a secular time scale in the presence of viscosity. The axisymmetry 
of the original figure is preserved but the topology (i.e. the connectedness of volume) is 
broken. The two tori along the lower branch of the two-ring sequence have inner-to-outer 
mass ratio of about 0.9. Other systems of two tori with different mass ratios do not bifurcate 
from the Maclaurin sequence; they seem to connect to various core-ring sequences (Eriguchi 
& Hachisu 1983a). 

Neither the bifurcation point nor the A-point of the two-ring sequence (Table 1) is 
related to the axisymmetric dynamical instability found by Bardeen (1971) at e = 0.998917, 
il = 4.598734, u? Q = 3.24702, E = -10.53929, and T/\W\ = 0.457424. This reinforces our 
conclusion in Paper II that Bardeen's dynamical instability represents a type 1 A-transition 
from the bifurcation point of the Maclaurin toroid sequence to the stable lower energy branch 
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of the same sequence. 

In its simplest level, the type 4 A-transition of the two-ring sequence is part of the "cusp 
catastrophe." Prior to the intersection X=C in Figure 3, the Maclaurin line represents stable 
equilibria while the two-ring branches represent maxima or inflection points of higher free 
energies. Such maxima and inflection points can only be explained at a higher level as parts 
of a higher elementary catastrophe. However, near point X=C, the free-energy diagram is 
simplified considerably. The two sequences exhibit degeneracy and their roles are reversed 
past this point. It is this smooth exchange of roles across X=C that can be described by 
an evolutionary path in the control plane of the cusp catastrophe. This path lies outside 
the fold lines that characterize doubly degenerate critical points (see Figure 5.4 and §5.4 
in Gilmore 1981). Staying outside the fold lines eliminates the first-order phase transition 
encountered prior to type 1 and type 2 A-transitions (Paper II). 

In summary, the above type 4 A-transition differs from the type 1 A-transition of Paper 
II in four respects: (a) it does not exhibit a first-order phase transition prior to the A-point; 
(b) it represents a secular instability; (c) the "specific heat" diverges at the A-point; and 
(d) it results in "double fission" to two toroidal figures. On the other hand, both types of 
transitions share several common properties since axisymmetry is preserved, topology breaks, 
and both evolutionary paths are parts of the cusp catastrophe (see also Table 3 below). 

4 Type 3 A-transitions 

Third-harmonic dynamical instability sets in at e=0. 96696 on the Maclaurin sequence 
(EFE) but there is no known sequence that bifurcates at that point. As we have briefly 
discussed in Paper II, the point of third-harmonic dynamical instability represents a type 
3 A-transition between the Maclaurin sequence and the lower energy stable branch of the 
dumbbell-binary sequence. This A-transition is shown in Figures 4 and 5 along with two 
more A-transitions of the same type that originate on the lower branch of the one-ring 
sequence. 

We have determined the properties of the type 3 A-points based on the j -values of the 
turning or starting points of the corresponding sequences. Analytical values of the relevant 
critical points on the Maclaurin sequence were obtained from Chandrasekhar (1968) and 
from EFE. Our results are summarized also in Table 2: 

(a) The A-transition M— >■ D— » B in Figures 4 and 5 was determined from the value j% = 1.540 
at the higher turning point D of the dumbbell-binary sequence. The A-point M with the same 
j -value on the Maclaurin sequence is identical to the point of third-harmonic dynamical 
instability (M' in Table 2). The transition leads to a detached binary (point B in Table 2 
and in Figures 4 and 5). 
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FIGURE 4. The rotation frequency squared along several incompressible sequences is 
plotted as a function of the dimensionless angular momentun squared The one-ring 
sequence was obtained from Table I of Eriguchi k, Sugimoto (1981). The dumbbell-binary 
sequence was obtained from Table II of Eriguchi, Hachisu, & Sugimoto (1982) and from 
Table 1 of Hachisu & Eriguchi (1984a). All points (open circles) of the triple and quadruple 
sequences were obtained from Table I of Eriguchi & Hachisu (1983b). The Maclaurin and 
Jacobi sequences and the bifurcation points (squares) of the dumbbell-binary and one-ring 
sequences (see Paper II) were calculated analytically. Point M lies directly above the higher 
turning point D of the dumbbell sequence where j% =1.540. The point of third-harmonic 
dynamical instability on the Maclaurin sequence with j„=l. 53796 (corresponding to a merid- 
ional eccentricity of e=0. 96696; EFE) coincides with M. Points M 3 and M 4 , corresponding to 
e=0. 98126 and 0.98495, lie directly above points R 3 and S4 (and R±), respectively. They were 
determined using the values j^=l. 980 (lower turning point R 3 of the one-ring sequence) and 
jg=2.159 (starting point S 4 of the quadruple sequence), respectively. Two of the three known 
fourth-harmonic neutral points (corresponding to e=0. 98097 and 0.98531; Chandrasekhar 
1968, EFE) essentially coincide with points M 3 and M 4 . Three type 3 X-transitions, leading 
to various multi-fluid-body sequences, are indicated by small solid dots (see also Table 2). 
Letter symbols indicate points discussed in the text. 
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FIGURE 5. The free energy E is plotted versus j% for the sequences of Figure 4. Three 
type 3 X-transitions, leading to various multi-ttuid-body sequences, are indicated by small 
solid dots (see also Table 2). Some letter symbols and labels and the open circles on the 
quadruple sequence seen in Figure 4 are omitted here for clarity. 

(b) The A-transition R 3 — >■ S 3 — > T in Figures 4 and 5 was determined from the value j% = 
1.980 at the lower turning point R3 of the one-ring sequence. An alternative determination, 
based on the value j 2 = 1.954 of the starting point S3 of the triple sequence, yields essentially 
the same answer. The transition leads to a three-fluid-body system (point T in Table 2 and 
in Figures 4 and 5). On the Maclaurin sequence, the above two j -values correspond to 
points M 3 and S' 3 , respectively, that are also listed in Table 2. These points are essentially 
both identical to the only fourth-harmonic neutral point (M 3 in Table 2) for which Eriguchi 
& Hachisu (1982) did not find a bifurcating sequence. 

(c) The A-transition R 4 — > S4 — >■ Q in Figures 4 and 5 was determined from the value 
jg = 2.159 at the starting point S4 of the quadruple sequence. This transition leads to a 
four-fluid-body system (point Q in Table 2 and in Figures 4 and 5). On the Maclaurin 
sequence, the above j -value corresponds to point M 4 that is also listed in Table 2. This 
point is esentially identical to the fourth-harmonic neutral point (M 4 in Table 2) which is 
also the bifurcation point of the one-ring sequence (see Table 1). 



Table 3 

Comparison Between A-Transitions^) 



Type 


Instability 


Spatial 
Symmetry 


Specific 
Heat 


First-Order 
Transition? 


Elementary 
Catastrophe 


Related 
Sequence 


Thermodynamics 
Analog^ 


1 


Dynamical 


Preserved 


Continuous 


Yes 


Cusp 


Maclaurin 
Toroid 


Bose-Einstein 
Condensation 


2 


Secular 


Preserved 


oo 


Yes 


(none) 


One-ring 
& Binary 


Superffuid 
Liquid 4 He 


3 


Dynamical 


Broken 


oo 


Yes 


7(c) 


Multi-fluid 
Bodies 


Binary Alloys 


4 


Secular 


Preserved 


oo 


No 


Cusp 


Two-ring 


Ferromagnetism^ 



(a) Topology breaks and fission occurs in all cases; see also Paper II. 

(6) See Huang (1963) for descriptions of the thermodynamical transitions. 

(c) If elementary, perhaps then symmetry-restricted butterfly catastrophe; see §9.5 in Poston & Stewart (1978). 

(d) See §14.2 in Poston & Stewart (1978). 
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All the above type 3 A-transitions represent dynamical instabilities that lead to break- 
ings of both the symmetry and the topology of the original figure. Related examples from 
linear stability analyses and from hydrodynamical simulations are discussed in §6 below. 
The breaking of the symmetry should take place prior to the breaking of the topology. This 
implies that type 3 A-transitions cannot take place before the corresponding critical points 
are reached because presumably there exist some sequences with branches of unstable equi- 
libria acting as free-energy barriers. Such unstable equilibria do not belong to the unstable 
branches of the multi-fluid-body sequences which have lower free energies than the mother 
(Maclaurin and one-ring) sequences (Table 2 and Figure 5). Nonetheless, the coincidence 
between critical points seen in Figures 4 and 5 suggests some kind of alignment between turn- 
ing or critical points in the two branches (the multi-fluid-body and the unknown branch) 
that represent unstable equilibria. 

It is not known which sequences act as free-energy barriers of type 3 A-transitions and, 
consequently, whether the barriers disappear smoothly or discontinuously at the A-points in 
Figures 4 and 5. Thus, it is not known whether type 3 dynamical transitions can be described 
by one of Thorn's (1975) elementary catastrophes. If the barriers disappear smoothly, then 
the so-called "symmetry-restricted butterfly catastrophe" (Poston & Stewart 1978, §9.5; 
Gilmore 1981, §10.7) is a good candidate for the description of these transitions because its 
control plane exhibits three distinct regions where the properties of the free-energy function 
differ. A system evolving along a path that crosses all three regions may then undergo 
symmetry breaking in the second region and topology breaking in the third region. 

A detailed comparison between the various types of hydrodynamical A-transitions is 
provided in Table 3 where we have collected the results of Paper II and of §3 and §4 of 
this paper. Table 3 summarizes our classification of A-transitions and relates them to the 
equilibrium sequences analyzed in this series of papers as well as to the corresponding ther- 
modynamical phase transitions and, where possible, to elementary catastrophes. The most 
interesting result contained in this table is that type 2 A-transitions are not elementary 
catastrophes because the free-energy barriers do not disappear smoothly at the correspond- 
ing A-points, suggesting the absence of degeneracy and the existence of a discontinuity at 
the critical points of the free-energy functions themselves (see also Paper II). 

5 Nonlinear resonances 

Figure 5 above was meant to illustrate the relative free energies between equilibria at 
the two ends of type 3 A-transitions. However, another interesting effect can also be seen 
in Figure 5. The lower turning point of the dumbbell-binary sequence and the turning 
points of the triple and quadruple sequences approach very close to the Maclaurin sequence 
in the (j%,E ) plane although they remain at slightly lower free energies relative to the 
Maclaurin points of the same angular momentum. This kind of structure in the free energies 
suggests the possibility of nonlinear resonances. The three nonlinear resonances of Figure 
5 are illustrated also in the (j%,u)o) plane of Figure 6 by small open circles. The physical 
properties of each resonance are summarized in Table 4. 
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FIGURE 6. The rotation frequency squared along the incompressible sequences of Figure 
4 is plotted again as a function of the dimensionless angular momentun squared j% but three 
nonlinear resonances between the Maclaurin sequence and the multi-fluid-body sequences 
are now indicated by small open circles. The resonances occur at the values j 2 =1 . 1 664, 1 . 704, 
and 2.082 of the turning points A, P, and C, respectively, of the multi-fluid-body sequences 
(see also Table 4). These values were obtained from the corresponding tables published by 
Eriguchi & Hachisu (1983b) and by Hachisu & Eriguchi (1984a). On the Maclaurin sequence, 
the resonances appear at points X, Y, and Z that correspond to values of the meridional 
eccentricity e=0. 94455, 0.97343, and 0.98348, respectively. The equilibrium free energies E 
in these resonances may be seen in Figure 5 which shows the turning points of the binary, 
triple, and quadruple sequences approaching very close to the Maclaurin sequence in the 
(jg, E Q ) plane. Letter symbols indicate points discussed in the text. 

The Maclaurin resonant points X, Y, and Z in Figure 6 and in Table 4 have only slightly 
higher free energies than the corresponding turning points of the multi-fluid-body sequences. 
The transitions to the lower energy states are not allowed to occur automatically because 
each pair of equilibrium states is presumably separated by a free-energy barrier. Such a 
free-energy barrier must be overcome by finite-amplitude, nonlinear perturbations before a 
Maclaurin spheroid may break in multiple bodies. So, the resonant points X, Y, and Z in 
Figure 6 may be thought as critical points of "first-order phase transitions." However, it is 
not known which equilibrium sequences play the role of the free-energy barriers. 
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Point X occurs at j% = 1.1664, i.e., prior to the point of second-harmonic dynamical 
instability on the Maclaurin sequence {j 2 = 1.27957; Table 1). Technically then, finite- 
amplitude disturbances can cause fission of the Maclaurin spheroid at X although it is dy- 
namically stable to all infinitesimal perturbations. On the other hand, points Y and Z occur 
past point M in Figure 6 when the Maclaurin spheroid is already linearly unstable to both 
second-harmonic and third-harmonic perturbations. A point of fourth-harmonic dynamical 
instability on the Maclaurin sequence has not been located yet, so we do not know whether 
it occurs prior to or past point Y or Z in Figure 6. 

We see from Figures 5 and 6 that the X — > A resonance leads to the stable branch of 
the dumbbell-binary sequence, the Y — > P resonance leads to the stable branch of the triple 
sequence, and the Z — > C resonance leads to the stable branch of the quadruple sequence. 
Let us now imagine nonlinear evolution driven by finite-amplitude disturbances and taking 
place at these resonances. Once at the lower turning point of the dumbbell-binary sequence, 
a system does not necessarily continue to evolve along the binary branch. It may just 
as well proceed to an even lower energy state on the Jacobi sequence by coalescence and 
further release of energy (see Figure 5). Similarly, once at the turning points of the triple 
or the quadruple sequence, a system may eject one or more of its components and make the 
transition to a state with fewer fluid bodies and lower energy. The ejection mechanism may 
break the entire multiple system apart to single unbound components but, in some cases, 
it may also lead to the formation of a bound binary system. In all cases, further evolution 
beyond the turning points becomes possible because of the strong nonlinear disturbances 
that have to be present in order to initiate this type of transition at a resonance. 

6 Discussion and Summary 

We have described the phase transitions that appear at higher than second-harmonic 
critical points on the Maclaurin sequence. These critical points are listed in Table 1. We 
have found that second-order phase transitions at the third-harmonic bifurcation of the tri- 
angle sequence and at the fourth-harmonic bifurcation of the square sequence are forbidden 
because the daughter sequences exist at higher energies relative to the Maclaurin sequence. 
Similarly, third-order phase transitions at the third-harmonic bifurcation of the ammonite 
sequence and at the axisymmetric sixth-harmonic bifurcation of the two-ring sequence are 
also energetically forbidden. In addition, the ammonite sequence exists only at higher rota- 
tion frequencies relative to the Maclaurin sequence and exhibits exactly the same structure 
as the pear-shaped sequence that bifurcates from the Jacobi sequence (Paper III). The trian- 
gle, square, and ammonite sequences terminate prematurely due to equatorial mass shedding 
while they are still unstable relative to the Maclaurin sequence (see Figure 1 in §2). Since 
there exist no lower energy branches, no instability of any kind (including A-transitions) 
is associated with the triangle, square, and ammonite sequences. These three sequences of 
objects in uniform rotation without any allowed phase transition are counter-examples to 
the common notion that bifurcation points are unconditionally points of secular instability 
as well. (See also Papers II and III for more related examples and for a discussion of this 
point.) 
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The two-ring sequence described in §3 has four turning points that change its structure 
and dynamical properties. Like the one-ring sequence in Paper II, the two-ring sequence 
does have a lower energy branch which however is accessible only via a new type of A- 
transition. We call this topology-breaking, symmetry-preserving transition a "A-transition 
of type 4" because, in contrast to the other types, it is not preceded by a first-order phase 
transition (see Figures 2 and 3 in §3). The type 4 A-transition of the two-ring sequence 
leads on a secular time scale to nonlinear fission to two coaxial axisymmetric toroidal figures 
that rotate uniformly and with the same frequency. 

We have also described two more kinds of nonlinear fission processes, the A-transitions of 
type 3 (§4, Table 2, and Figures 4, 5) and the nonlinear resonances of the Maclaurin sequence 
(§5, Table 4, and Figures 5, 6). The point of third-harmonic dynamical instability on the 
Maclaurin sequence lies exactly at the same angular momentum value as the higher turning 
point of the dumbbell-binary sequence. Thus, this point appears to be a A-point. The A- 
transition is illustrated in Figures 4 and 5. It is of type 3 because the symmetry is also broken 
in addition to the topology. In fact, the symmetry breaks first since the linear dynamical 
instability is caused by third-harmonic perturbations. The topology breaks second during 
nonlinear evolution and this A-transition finally results in fission of a Maclaurin spheroid to 
a binary. This type of evolution has been discussed in Paper II and has been observed in the 
numerical, smoothed-particle hydrodynamical simulations of Miyama, Hayashi, & Narita 
(1984) during nonequilibrium evolution of a hot isothermal model with T/|W| = 0.3. This 
value is in good agreement with the critical value of T/|W| = 0.30329 listed in Table 2. 

We have presented in §4 two additional examples of type 3 A-transitions that originate 
on the lower toroidal branch of the one-ring sequence. The A-points on the one-ring sequence 
appear to be related to the starting points of the triple and quadruple sequences, respectively, 
as well as to two of the three known neutral fourth-harmonic points on the Maclaurin 
sequence. These relationships are illustrated in Figures 4 and 5 and are quantified in Table 
2. It is not known which equilibrium sequences act as free-energy barriers preventing the 
transitions prior to the A-points of type 3. 

The physical characteristics of the four types of A-transitions are compared in Table 3. 
The analogies drawn between the hydrodynamical A-transitions and the thermodynamical 
phase transitions and catastrophes are based on the common properties described in §2 and 
§3 above and in Paper II. However, these analogies are speculative for the most part and 
need to be confirmed by future work. 

The type 3 A-transitions from the one-ring sequence lead on a dynamical time scale 
to fission to three and four detached fluid bodies. These dynamical instabilities are impor- 
tant because their exact counterparts have been discovered by linear stability analyses of 
(in)-compressible, differentially rotating, self-gravitating tori and annuli (this confirms once 
again the detection of dynamical type 3 A-transitions by linear analyses) and have been 
thoroughly studied in the nonlinear regime with the help of numerical simulations. Such 
unstable modes were first discovered by Goodman & Narayan (1988) who called them the 
"Intermediate" or "I-modes." Subsequent linear analyses and numerical simulations of two- 
dimensional annuli helped to identify the I-modes with nonlinear fission (Christodoulou &: 
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Narayan 1992; Christodoulou 1993). The nonlinear manifestations of the instability were also 
studied numerically in three-dimensional tori (Tohline & Hachisu 1990; Woodward, Tohline, 
& Hachisu 1994). Our results leave no doubt that the dynamically unstable I-modes in tori 
and annuli are responsible for nonlinear fission to multiple fluid bodies and, in addition, 
they clearly identify these I-modes with type 3 A-transitions. Combined with the analyt- 
ical and the numerical work on I-modes cited above, this identification suggests that the 
A-transitions discussed in this series of papers do occur also in compressible, differentially 
rotating, self-gravitating fluids and are thus relevant to the evolution of astrophysical fluid 
systems as well. 

We have demonstrated in §5 the power of the nonlinear method that we employ to study 
instabilities by using the equilibrium free energies shown in Figure 5 to identify nonlinear 
resonances between the Maclaurin sequence and various multi-fluid-body sequences. Reso- 
nances occur at points where the free energies of the binary, triple, and quadruple sequences 
are only slightly lower than those of the Maclaurin sequence for the same angular momentum 
value. We have described the physical properties of three such resonances based on Figures 
5 and 6 and we have summarized the relevant parameters in Table 4. Resonances are fully 
nonlinear phenomena and cannot be detected by stability analyses that are accurate to any 
particular order. Although the free-energy differences are small at resonances (see Figure 5) , 
phase transitions such as those depicted in Figure 6 can only be induced by finite-amplitude 
disturbances because free-energy barriers presumably stand between the two equilibrium 
states. Therefore, resonances effectively appear as critical points of "first-order phase tran- 
sitions" since a lower energy state has opened up to an evolving Maclaurin spheroid but a 
free-energy barrier prevents the transition from taking place automatically. It is not known 
which equilibrium sequences form the free-energy barriers. 

Of the three nonlinear resonances listed in Table 4 and illustrated in Figure 6, the 
Maclaurin-to-binary resonance is the most interesting because it appears at T/|W| = 
0.25947, i.e. prior to the point T/|W| = 0.27382 of second-harmonic dynamical insta- 
bility. This last point is the lowest point of dynamical instability on the Maclaurin sequence 
predicted by linear analysis. Thus, a Maclaurin spheroid subject to finite-amplitude distur- 
bances at this resonance may break up to a binary although it is dynamically stable to all 
infinitesimal perturbations. 

As we have mentioned in §5 above, because of the extreme nonlinearities involved, 
dynamical evolution at each resonance does not necessarily have to continue along the binary, 
triple, or quadruple sequence, respectively, since there exist nearby states of even lower 
energy. Specifically, a multiple fluid system produced by finite-amplitude disturbances may 
coalesce forming a Jacobi ellipsoid or it may even eject some/all of its components producing 
a lower energy system with fewer components or several single unbound bodies (cf. Figures 
5 and 6). In some cases, the ejection mechanism may however produce a bound binary, so 
this process may be of some importance to the formation of binary stars. 

Equilibrium sequences of multiple fluid bodies with different compressibilities as well as 
binary sequences with different mass ratios have been calculated numerically by Eriguchi & 
Hachisu (1983b, 1984), Hachisu & Eriguchi (1984a, b), and Hachisu (1986) who also discussed 
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briefly some of the phase transitions between sequences. Their results help us understand 
better the nonlinear fission processes that we have described mainly in this paper and in 
Papers II and III. Such processes may be responsible for the formation of multiple stellar 
systems and, in particular, binary stars that seem to be more common than single stars in 
our Galaxy. [See the reviews of Abt (1978, 1983) and Bodenheimer (1992), the recent work 
of Abt, Gomez, & Levy (1990), and references therein.] Capture of stars in stellar clusters 
has been proposed as an alternative mechanism of binary-star formation. It is, however, 
unlikely that capture may be responsible for the formation of the observed short-period, 
spectroscopic binaries. 

The capture mechanism strives to explain the observed unequal-mass binaries (Abt, 
Gomez, & Levy 1990). The presumption is that the classical fission mechanism from the 
Jacobi sequence (Paper III) can only lead to the formation of equal-mass binaries. Although 
this is true, capture is not the only mechanism that produces unequal-mass binaries. The 
evolutionary paths described above, and the type 3 A-transitions in particular, may also 
end to unequal-mass binaries. This is particularly easy to achieve via the A-transitions 
that originate from the toroidal states of the one-ring sequence. These transitions may 
easily form unequal-mass multiple systems especially when the original toroidal structures 
are unstable to several modes of comparable growth rates (cf. Christodoulou & Narayan 
1992). Unequal-mass binaries may then be left behind after the ejection of the additional 
low-mass components. Therefore, the presumption that unequal-mass binaries can form 
only by capture is no longer justified in view of such nonlinear fission processes. 

But how can unequal-mass short-period binaries form? It is difficult to have such 
systems form by capture and it is not clear whether the above nonlinear fission processes can 
lead to their formation. The unequal-mass binary sequences lie between the stable parts 
of the Maclaurin and Jacobi sequences and the equal-mass binary sequence (Hachisu & 
Eriguchi 1984b; see also Paper III). Perhaps these sequences can be accessed from the stable 
parts of the Maclaurin and the Jacobi sequences through additional nonlinear resonances 
but this matter has not yet been explored. 

At this point, we close our discussions of results as well as this entire series of papers with 
an important reminder: Most of our results and conclusions have been obtained by analyzing 
rotating, self-gravitating, incompressible fluid models and cannot be extended by default to 
the studies of compressible and/or magnetized systems where internal, magnetic, and other 
forms of energy may play an important role. (Such energy sources have been excluded from 
the free-energy functions that we have used so far.) The theory of phase transitions is, 
however, applicable also to compressible (see Paper III) and to stellar (see Paper I) systems 
and is expected to yield a clearer physical interpretation of the instabilities, catastrophes, 
and bifurcations encountered in such systems as well. 
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Table 2 
Type 3 A-Transitions 
Point jl g E T/\W\ e 

MACLAURIN TO BINARY SYSTEM: M -> D ->• B 

~~ M 1.540 10.48661 -7.10825 0.30329 0.96705 
D 1.540 3.288 -7.754 0.18260 
B 1.540 0.3594 -8.002 0.06700 
Third-harmonic Dynamical Instability (EFE) 
"M* 1.53796 10.49141 -7.10294 0.30308 0.96696 
ONE-RING TO TRIPLE SYSTEM: R 3 — S 3 — )► T 



R 3 


1.980 


7.556 


-8.031 


0.32280 


S 3 


1.954 


3.455 


-8.430 


0.22900 


T 


1.980 


1.044 


-8.612 


0.14230 



On the Maclaurin Sequence 

~M 3 1.980 9.29812 -8.10289 0.34393 0.98126 

S 3 1.954 9.37373 -8.05183 0.34180 0.98065 
Fourth-harmonic Neutral Point (EFE) 

~m£ 1.96741 9.33478 -8.07828 0.34290 0.98097 

ONE-RING TO QUADRUPLE SYSTEM: R 4 S 4 ->• Q 



R4 


2.159 


5.796 


-8.417 


0.31200 


s 4 


2.159 


3.826 


-8.490 


0.26800 


Q 


2.159 


2.324 


-8.514 


0.21960 



On the Maclaurin Sequence 



M 4 2.159 8.77090 -8.43155 0.35780 0.98495 
Fourth-harmonic Neutral Point (EFE) 



M 4 2.17884 8.71210 -8.46563 0.35925 0.98531 
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Table 4 
Nonlinear Resonances 

Point f g E T/\W\ e 

MACLAURIN TO BINARY SYSTEM: X -)• A 

1.1664 11.15073 -6.00275 0.25947 0.94455 
A 1.1664 1.765 -6.175 0.11800 
MACLAURIN TO TRIPLE SYSTEM: Y -)• P 

1.704 10.07324 -7.51326 0.31964 0.97343 
P 1.704 2.853 -7.620 0.19800 
MACLAURIN TO QUADRUPLE SYSTEM: Z -)• C 



Z 2.082 8.99873 -8.29493 0.35200 0.98348 
C 2.082 3.323 -8.300 0.24800 



